Scaling laws in the quantum to classical transition in chaotic systems 
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We study the quantum to classical transition in a chaotic system surrounded by a diffusive en- 
vironment. The emergence of classicality is monitored by the Renyi entropy, a measure of the 
entanglement of a system with its environment. We show that the Renyi entropy has a transition 
from quantum to classical behavior that scales with h^g/D, where h e a is the effective Planck con- 
stant and D is the strength of the noise. However, it was recently shown that a different scaling law 
controls the quantum to classical transition when it is measured comparing the corresponding phase 
space distributions. We discuss here the meaning of both scalings in the precise definition of a fron- 
tier between the classical and quantum behavior. We also show that there are quantum coherences 
that the Renyi entropy is unable to detect which questions its use in the studies of decoherence. 



Since the birth of the quantum theory, the fundamen- 
tal question of how a quantum mechanical system starts 
to behave classically has been the subject of an intense 
debate. In the last years, it has been well established 
that decoherence induced by the environment is the main 
ingredient in the quantum-to-classical transition [l|, Q. 
However, a clear way to describe the fuzzy boundary be- 
tween the quantum and the classical world has not been 
found yet. 

Many attempts have been done to understand the loss 
of the quantum coherence and several measures of the 
restoration of the classicality have been proposed spe- 
cially in chaotic systems. Recently, it was conjectured 
that, in the presence of noise, a single composite param- 
eter controls the quantum-to-classical transition of clas- 
sically chaotic systems This single parameter should 
appear in the computation of some measure which di- 
rectly reflects "the distance" between quantum and clas- 
sical evolution as a function of the effective Planck con- 
stant K e ff {i.e the relative size of the Planck constant), 
the strength D of the coupling with the environment and 
the Lyapunov coefficient A of the classical dynamics of 
the system. Thus, the computed "distance" would be 
controlled by a combined parameter of the general form 
£ = Ti" s \PD-f when the transition from quantum to clas- 
sical behavior occurs. This conjecture was shown to be 
valid, for example in 0] , where the phase space integral of 
the modulus of the difference between the Wigner func- 
tion and the classical distribution was used as a measure 
of the quantum to classical "distance" . 

Another usual quantity to monitor the transition from 
quantum to classical behavior is the purity P = Tr[p 2 ], 
which is a measure of the entanglement of the sys- 
tem with the environment. Any function of the purity 
has the same information but the so called Renyi en- 



tropy, S = — log(P), is particularly interesting because 
it clearly reveals that above a threshold of the strength of 
the coupling with the environment its behavior is dom- 
inated by the positive classical Lyapunov exponents of 
the system Q . 

In this work we show, for the same system and environ- 
ment used in [ij , that the behavior of the Renyi entropy 
S (and thus the purity P) is controlled by the combined 
parameter \' = fi-efi/D- This scaling was also founded 
in 0] using the Duffing oscillator, suggesting an univer- 
sality of this scaling law. In our system the emergence 
of classicality when monitored with the Renyi entropy 
should correspond to the regime \' < 100. However, 
in the regime 1 < "x! < 100 the evolved Wigner func- 
tion displays an interference fringe pattern that clearly 
indicates that not all the quantum coherences of the sys- 
tem have dissapeared. We explain that while \' controls 
which quantum terms effectively contribute in the mas- 
ter equation for evolution of the Wigner function, the 
fine coherent structure of this function is controlled by 
the composite parameter \ = We also 

show that only in this essentially classical regime and for 
short times, the behavior of the Renyi entropy is inde- 
pendent on the diffusion constant and dominated, in our 
case, by a classical Lyapunov exponent. 

The system we analyze is the kicked harmonic oscilla- 
tor (KHO) whose quantum Hamiltonian is 

P 2 1 

H = h -mv 2 Q 2 + A cos(kQ) d n (r), (1) 

2m 2 

where 5 n {r) = X)^o — nT )^ m ^ s the particle mass, 
v the oscillator frequency, A the amplitude of the kicks 
and r the interval between two consecutive kicks. Work- 
ing with dimcnsionlcss quantities q = kQ, p = kP/mv 
and K = k 2 A/mv, we can define an effective Planck con- 
stant h e s, so that [q,p] = 2ir] 2 = ifi c s, with the so called 
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FIG. 1: Quantum (•) and classical (o) Renyi entropy as a 
function of the nu mber of the kicks for the KHO with K = 2 
and rj = \Jh e g/2 — 0.125. From bottom to top, the diffu- 
sion constants are D = 2.56 x 1CT 7 ,6.4 x 1CT 6 ,9 x 1CT 5 and 
4 x 10 -3 . The dashed line is a linear function with a slope 
given by the logarithm of the classical expansion coefficient 
at (q,p) = (0,0). 



Lamb-Dicke parameter rj = kAQo = k^fi/2mv being 
the quotient between the width AQq of the ground state 
of the harmonic oscillator and the wavelength A = 2-7r/fc 
of the position dependent pulse. The classical dynamics 
in phase space of the KHO is unbounded and the stro- 
boscope sections exhibit a mixed dynamics with stable 
islands surrounded by a chaotic sea Q. Here, we con- 
sider the relation r = T/6 between the period of the 
harmonic oscillator T = 2ixjv and the interval between 
kicks r. Thus, the classical phase space structure corre- 
sponds to a stochastic web, its thickness governed by K. 
We have considered K — 2. thus the phase space origin 
is an hyperbolic fixed point and the area around it shows 
an essentially chaotic behavior 

In order to study the decoherence effects that induce 
a quantum-to-classical transition we coupled the KHO, 
as usual 0, 0, @, , to a thermal reservoir with average 
population n and coupling constant T in the Markovian 
and weak coupling limit, and we considered the purely 
diffusive regime (i.e. the high temperature regime n — > 
oo together with T — > and T = nTrf a constant value). 
For this kind of environment the master equation for the 
evolution of the Wigner function W(x) reads 



dW 
~dt 



L cl + L q +T= {H, W} 
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The equivalent evolution equation for the classical distri- 
bution only contains the evolution generated by the Pois- 
son bracket L c i and the term T that comes from the cou- 
pling with the diffusive environment. The quantal term 
L q comes from the nonlinear part of the Hamiltonian that 
in our case are simple the kicks V(q) = K cos(g)5 n (r). 
The strength of the interaction with the enviroment is 
controled by the dimensionless diffusion constant, D = 
fr. 

Let us first consider the quantum-to-classical tran- 
sition by monitoring the behavior of the classical 
and quantum Rcnyi entropies S = — log(P) = 
— log(27r?i c ff J W 2 (x)dx) (W is the Wigner function or 
the classical phase space distribution). In Fig. Q]we show 
the behavior of the entropies in the KHO with K = 2 and 
7] = \jTi c {{/2 = 0.125 for several values of the strength of 
the interaction with the enviroment D. The initial state 
is a coherent state centered at the origin of phase space 
with width Aq = Ap = rj. We can see that for extremely 
small values of D (D = 2.56 x 10 -7 ), and short times, the 
quantum and classical evolution are nearly unitary, and 
thus the Renyi entropies are S n ~ 0. Increasing the value 
of the diffusion constant D and in the regime where the 
initially localized wave packet stretches in the direction 
of the unstable classical manifold the quantum and clas- 
sical Renyi entropy exhibit different time regimes [Fig. 
[p. When D is still relative small, we can see two regimes 
[l 21 ] . In the first one the behavior of both entropies is 
essentially quadratic in time and the differences between 
them are negligible. Then, there is a linear regime whith 
a slope that depends on D and can be different for the 
quantum and the classical case. This linear regime oc- 
curs for times of the order of the well known log-time 
t& = log(l/7i c ff)/(2A), with A the Lyapunov exponent 
or, as in our example, the logarithm of the expansion co- 
efficient of the linearization of the dynamics at the origin 
(hyperbolic fixed point) 0]. Finally, we can see that for 
sufficiently greater values of D the linear growth of the 
quantum and the classical Renyi entropy are identical 
and given by A 0]. 

We have studied systematically the dependence of the 
linear growth regime with the diffusion constant D and 
h c g. In Fig. (inset) we show the slope 7 for the 
linear growth of the quantum and classical Renyi en- 
tropy as a function of the diffusion constant D for sev- 
eral values of ?i c ff . Both slopes 7 approach A (horizontal 
line) for all the h e s considered but at different values 
of the diffusion constant D. However, when 7 is plot- 
ted as a function of h eS /D a clear scaling emerges. We 
can see that the classical and quantum entropy linear 
growth rate measured by 7 are equal for h eS /D < 100, 
but it is only in the range fi 2 s /D < 1 that the classi- 
cal and quantum slope approach A (the so called "Lya- 
punov" regime). On the other hand, for h eS /D > 100 
the linear growth rate of the Renyi entropy depends 
linearly with the diffusion constant D. This is the so 
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FIG. 2: Slope 7 of the linear growth of the quantum (full 
symbols) and classical (empty symbols) Renyi entropy as a 
function of h^/D (main plot) and 1/D (inset) (D is the dif- 
fusion constant) for 77 = yJh c s/2 = 0.5 (triangles), 0.3125 
(diamonds), 0.125 (squares) and 0.0625 (circles). We have 
also plotted the functions 7 oc D/h^g (dotted line) and 
7 oc \ /r D/h c g (dashed line). See text for details. 
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FIG. 3: The measure a(riE) (see text) of the differences be- 
tween the quantum and classical Renyi entropy. The initial 
state is a coherent state centered around the origin of phase 
space with width Aq — Ap = r\, with r\ = v2^eff = 0.5 (trian- 
gles), 0.3125 (diamonts), 0.125 (squares) and 0.0625 (circles). 
For each h e ff the Renyi entropies were computed at the log- 
time £e w uet. 



called Fermi golden-rule regime that was explained in 
Refs. @, E3. For the classical case, we found that 
the Renyi entropy linear growth rate scales with \f~D. 
This fact deserves further investigation [l2| . The scaling 
with Ti 2 cfi /D can also be appreciated in Fig.© where we 
plot the measure, <j{n) = \S q (n) — S c i(n)\/S q (n), of the 
"distance" between the quantum Renyi entropy, S q , and 
the classical one, S c i, computed at the number of kicks 
n = n E w t E /r = log(l/?i eff )/(2rA). 

Another approach to monitor the emergency of classi- 



cality is to directly compare the evolution of the Wigner 
function and the corresponding classical distribution in 
the presence of the environment 0, HI ■ A different scaling 
law appears in this approach where it was shown that the 
Wigner function and the classical distribution are essen- 
tially equal in the parameters regime h 2 s /D << 1 and 
X = h 2 s K/4:D 3 / 2 < 1. However, we have seen that the 
linear growth of the classical and quantum Renyi entropy 
are equal for h^g/D < 100. Even more, from Fig. ([I]) we 
see that the differences between the quantum and classi- 
cal Renyi entropy are negligible beyond the linear regime 
up to n < 9 for example for D = 9 x 10~ 5 where x' = 
hl s /D ks 10 < 100 but x = hl^K/AD 3 / 2 w 571 > 1. In 
Fig. [U we show the quantum and classical distribution 
immediately before the kick n = 7 for x' ~ 10 where we 
clearly see that the Wigner distribution shows important 
interference fringes arising from quantum coherences not 
detected by the Renyi entropy. These quantum coher- 
ences separate the quantum from the classical behavior. 
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FIG. 4: (Color online) Wigner distribution (top panel) and 
classical distribution [inset (a)] inmediatly before the kick n = 
7 for 77 = \/2h c g — 0.125 and difussion constant D = 9 x 10 . 
Inset (b) Wigner distribution inmediatly before the kick n = 8 
for 77 = \/2h f[ = 0.0256 and difussion constant D = 9 x 10~ 5 . 
Bottom panel: Wigner distribution (blue dotted line) and 
classical distribution (red solid line) evaluated on (q=-2,p) 
(dashed lines in the top panel). 

We can understand the different role of the scal- 
ings x' and x from the exact solution of Eq. ([2]) be- 
tween kicks. This solution can be written as the map: 
W„+i(x) = / dx' L(x fl ,x') W„(x') where W n+ i and W n 
are the Wigner functions or the classical phase space dis- 
tributions immediately before the kicks n + 1 and n re- 
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spectively, L is the quantum or classical propagator of 
one kick plus the action of the diffusion reservoir and 
x i?. = ^qR„R\ = R _1 (x) correspond to the harmonic 
evolution [J, |9|. It is instructive to write the one step 
quantum nonunitary propagator L for an arbitrary posi- 
tion dependent kick V(q): 



L(x fl ,x') 



e 



dfx 



X 



VinD J-oo 2nh eS 

e - K ^[V{q'-^/2)-V{ q '+ fl /2)-^p R -p')\ ,^ 



(for V(q) = K cos(q) this formula corrects for a small 
mistake in Eq. (16) of Ref. 0] ). If we write, 



V(q'_)-V(q + ) = -fxd q V- £ 



,,2m+l 

J± d 2m+1 V 

- n (2m + l)\ q 

rn—l K 7 

(4) 

where q'± = </±u/2 and d 2m+1 V = d 2m+1 V {q ) / dq 2m+1 
(to = 0,1,...) |l3[, and we kept only the first term 
we recover the classical propagator of the classical dis- 
tribution: L cl (x R ,x') = e~( x2+y2 } /4nD where y 

V + d q V^j /2VD and x = (q R - q) /2y/D. Thus, 

we see from Eq.Q that the same high order odd deriva- 
tives that enter in the quantal "correction" L q in the 
master equation ([2]) enter in the one step propagator ((3}. 
Therefore, the value of the variance x'/2 = hl s /2D of 
the Gaussian factor exp(— D^ 2 /fi 2 ^ ) in Eq. (J3j) is cru- 
cial in the quantum-to-classical transition because it 
controls which terms of this quantal "correction" effec- 
tively contribute. Thus, if the variance x'/2 is suf- 
ficiently small for the quantal "corrections" we have 
\^\ 2m+1 \d 2m+1 V\/h cS c < \n\ 3 \dV max \/h cS c where c = 
2 2m (2m + 1)! and \d 2m+1 V\ < \dV max \ is an upper 
bound valid in the region where the initial wave packet 
spread (for the KHO, I^VmaxI = K in all position space). 
But, |/^| 3 |9Vi nax |/?i c ffC it is only a small correction iff 
X = h eS I^Vmaxl/^D 3 / 2 << 1 where, following the steps 
described in [9( , we can write 



L(x fl ,x') f«r'(x K ,x') 



i + d " V{q) ft ^ 

|ov; nax | 



(5) 



with f (y) = 1/4 (y - 2y 3 /3). Based on this result, 
in Ref. [1, 0] it is showed that the distance be- 
tween the corresponding phase-space distributions,D n = 
/ dx | W n (x) — W£ l (x) | , that quantify the quantum ef- 
fects, scales with the single parameter \. In Fig. Q 
(inset b) we can see that for x = 1 the evolved Wigncr 
distribution does not display essentially any quantum co- 
herence. 

We have shown that the quantum-to-classical transi- 
tion in the presence of a diffusive enviroment is governed 



by the single parameter x' = ft 2 /D when it is measured 
with the Rcnyi entropy. Also, we have shown that x' 
controls whose terms are the leading order quantum cor- 
rections in the master equation for thc Wigncr function. 
However, for a class of kicked systems we have shown that 
the combine parameter that indicates in which regime the 
leading correction is negligible is x = ^ c ff 

\dV max \/4D 3 / 2 

(that for the KHO reduce to x = ^IffK/^D 3 / 2 ). It was 
previously shown that the distance between the classical 
and quantum phase space distributions, which measure 
thc quantum signatures, is proportional to the parameter 
X = hl s K/4:D 3 / 2 . Thus, it is thc combined x parame- 
ter that defines the precise frontier between the classical 
and quantum behavior of a system in a purely diffusive 
reservoir. Furthermore, our results clearly indicate that 
the dynamics of the Renyi entropy or the purity is unable 
to detect not negligible quantum coherences and there- 
fore questioned its use to study the quantum to classical 
transition. 
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